Hydrodynamic Coupling of Two Brownian Spheres to a Planar Surface 
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We describe direct imaging measurements of the collective and relative diffusion of two colloidal 
spheres near a flat plate. The bounding surface modifies the spheres' dynamics, even at separations 
of tens of radii. This behavior is captured by a stokeslet analysis of fluid flow driven by the spheres' 
and wall's no-slip boundary conditions. In particular, this analysis reveals surprising asymmetry in 
the normal modes for pair diffusion near a flat surface. 



Despite considerable progress over the past two cen- 
turies |lj hydrodynamic properties of all but the simplest 
colloidal systems remain controversial or unexplained. 
For example, velocity fluctuations in sedimcnting col- 
loidal suspensions are predicted to diverge with system 
size Experimental observations indicate, on the other 
hand, that long-wavelength fluctuations are suppressed 
by an as- yet undiscovered mechanism ||-|| . One possible 
explanation is that hydrodynamic coupling to bounding 
surfaces may influence particles' motions to a greater ex- 
tent and over a longer range than previously suspected 
H . Such considerations invite a renewed examination of 
how hydrodynamic coupling to bounding surfaces influ- 
ences colloidal particles' dynamics. 

This Letter describes an experimental and theoretical 
investigation of two colloidal spheres' diffusion near a 
flat plate. Related studies have addressed the dynam- 
ics of two spheres far from bounding walls |?]||, and 
of a single sphere in the presence of one or two walls 
||. Confinement by two walls poses particular difficul- 
ties since available theoretical predictions apply only for 
highly symmetric arrangements fio)] , or else contradict 
each other Jll],[lJ] . The geometry we have chosen avoids 
some of this complexity while still highlighting the range 
of non-additive hydrodynamic coupling in a many-surface 
system. 

We combined optical tweezer manipulation Jl3| and 
digital video microscopy |H| to measure four components 
of the pair diffusion tensor for two colloidal spheres as a 
function of their center-to-center separation r and of their 
height h above a planar glass surface. Measurements 
were performed on silica sphere of radius 0.495±0.025 ^m 
(Duke Scientific lot 21024) dispersed in a layer of water 
140 ± 2 /im thick. The suspension was sandwiched be- 
tween a microscope slide and a #1 coverslip whose sur- 
faces were stringently cleaned before assembly |lq] and 
whose edges were hermetically sealed with a uv cured 
epoxy (Norland type 88) to prevent evaporation and sup- 
press bulk fluid flow. A transparent thin film heater 



bonded to the microscope slide and driven by a Lakeshore 
LM-330 temperature controller maintained the sample 
volume's temperature at T = 29.00 ± 0.01°C, as mea- 
sured by a platinum resistance thermometer. The ad- 
dition of 2 mM of NaCl to the solution minimized elec- 
trostatic interactions among the weakly charged spheres 
and glass surfaces by reducing the Debye screening 
length to 2 nm. Under these conditions, the individ- 
ual spheres' free self-diffusion coefficients are expected to 
be D = fc B T/(67r?/a) = 0.550 ± 0.028 /zmf/sec, where 
7/ = 0.817 cP is the electrolyte's viscosity 
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FIG. 1. Measurement of a typical pair diffusion coefficient 
in the geometry depicted in the lower inset. Dashed spheres 
represent hydrodynamic images. The upper inset shows the 
histogram of two spheres' collective displacements in the ± 
direction starting from h = 25.5 ±0.7 fim, r = 7.00 ±0.25 ^m, 
and p = 0, after free diffusion for r = 1/30 sec. Fitting to a 
Gaussian yields the rms displacement Ap±(r, h, r). The main 
plot tracks the evolution of Ap5_ (t, h, r) together with a least 
squares fit to Eq. ([j]) for the diffusion coefficient (r, h) . 
The result is indicated by an arrow in Fig. ^(d). 

The spheres' motions were tracked with an Olympus 
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IMT-2 optical microscope using a 100 x NA 1.4 oil immer- 
sion objective. Images acquired with an NEC TI-324A 
CCD camera were recorded on a JVC-822DXU SVHS 
video deck before being digitized with a Mutech MV- 
1350 frame grabber at 1/60 sec intervals. Field-accurate 
digitization was assured by interpreting the vertical in- 
terlace time code recorded onto each video field. The 
spheres' locations r\{t) and ^(i) in the image acquired 
at time t then were measured to within 20 nm using a 
computerized centroid tracking algorithm |3. 

A pair of spheres was placed reproducibly in a plane 
parallel to the glass surfaces using optical tweezers [ pi| |. 
These optical traps were created with a solid state laser 
(Coherent Verdi) whose beam was brought to a focus 
within the sample volume by the microscope's objec- 
tive. Resulting optical gradient forces suffice to localize 
a silica sphere at the focal point despite random thermal 
forces ]ll| . Two optical traps were created by alternating 
the focused laser spot between two positions in the focal 
plane at 200 Hz using a galvanometer-driven mirror [fl7| . 
Diverting the trapping laser onto a beam block every few 
cycles freed the spheres to diffuse away from this well 
defined initial condition. Resuming the trap's oscillation 
between the two trapping points resets the spheres' po- 
sitions. Alternately trapping and releasing the spheres 
allowed us to sample their dynamics efficiently in a par- 
ticular geometry. Allowing the spheres only r = 83 msec 
(5 video fields) of freedom before retrapping them for 
16 msec (less than 1 video field) ensured that their out-of- 
plane motions, Az < \/2Dqt = 0.4 /im, cause negligible 
tracking errors. 

Because optical tweezers form in the microscope's fo- 
cal plane, their height h relative to the coverslip's sur- 
face can be adjusted from 1 to 30 (im with 0.5 /nm ac- 
curacy by adjusting the microscope's focus. For a given 
height, we continuously varied the spheres' initial sepa- 
ration between 2 fim and 10 /im at 0.025 Hz for a total 
of 20 minutes. This procedure yielded 60,000 samples 
of the spheres' dynamics in 1/60 sec intervals divided 
into sequences 5/60 sec long for each value of h. These 
trajectory data were decomposed into cooperative mo- 
tions p = r\ + r-i and relative motions r = f\ — r*2 either 
perpendicular or parallel to the initial separation vector, 
and binned according to the initial separation, r. The 
diffusion coefficients (r, h) associated with each mode 
of motion ij)(r, h, r) at each height and initial separation 
were then obtained from the Stokes-Einstein formula 

(A^ 2 (r)) =2D^(r,h)T, (1) 

where the angle brackets indicate an ensemble average. 

Fig. |l| shows typical data for one mode of motion at 
one height and starting separation. Diffusion coefficients 
Dip (r, h) extracted from least squares fits to Eq. (Q) ap- 
pear in Fig. Has functions of r for the smallest and largest 
accessible values of h. The horizontal dashed lines in 
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FIG. 2. Pair diffusion coefficients for 1 /im diameter sil- 
ica spheres as a function of center-to-center separation r and 
at two different center-to-surface heights h. (a) and (b): 
h = 1.55 ± 0.66 /nm. (c) and (d): h = 25.5 ± 0.7 /xm. (a) 
and (c): Collective and relative motions parallel to the ini- 
tial separation vector, (b) and (d): Perpendicular. Dashed 
curves result from linear superposition of drag coefficients. 
Solid curves result from the theory described here, with no 
adjustable parameters. Horizontal dashed lines indicate the 
asymptotic self-diffusion coefficient D xy (h) from Eq. (0). 

Fig. H indicate the spheres' asymptotic self-diffusion co- 
efficients. Measurements' deviations from these limiting 
values reveal the influence of the spheres' hydrodynamic 
interactions. 

Particles moving through a fluid at low Reynolds num- 
ber excite large-scale flows through the no-slip bound- 
ary condition at their surfaces. These flows couple dis- 
tant particles' motions, so that each particle's dynam- 
ics depends on the particular configuration of the en- 
tire collection. This dependence is readily calculated us- 
ing Batchelor's generalization of Einstein's classic argu- 
ment : The probability to find N particles at equilib- 
rium in a particular configuration {ri, . . . , r/v} depends 
on their interaction $(ri, . . . , fjv) through Boltzmann's 
distribution, P(r*i, . . . , r/v) = exp [— $/(/cbT)]. The cor- 
responding force — V$ = ksTVP/P drives a proba- 
bility flux ksTh VP, where h(fi, . . . , rjv) is the parti- 
cles' mobility tensor. The system reaches equilibrium 
when this interaction-driven flux is balanced by a diffu- 
sive flux — DVP. It follows that the iV-particle diffusiv- 
ity is D = ksTh. Elements of D lead to generalized 
Stokes-Einstein relations {Tg[| 

(Ar« (t) Ar Jl3 (r)) = 2D la . jP r. (2) 

describing how particle i's motion in the a direction cou- 
ples to particle j's in the (3 direction. 
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The mobility tensor for spheres of radius a has the form 

&ict,j{3 



(3) 



bf„ a a is the Green's function for the flow at fi in the a 
direction due to an external force at rj in the beta di- 
rection. In the present discussion, it accounts for no-slip 
boundary conditions at all other surfaces in the system. 

If the spheres are well separated, we may approximate 
the flow field around a given sphere by a stokeslet, the 
flow due a point force at the sphere's location. This ap- 
proximation is valid to leading order in the spheres' ra- 
dius. The Green's function for the flow at x in the a 
direction due to a stokeslet at fj in the f3 direction is p3] 
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In the particular case of 



two identical spheres, diagonalizing the resulting diffu- 
sivity tensor D yields the diffusion coefficients for two 
collective (C) modes and two relative (R) modes along 
directions perpendicular (_L) and parallel (||) to the ini- 
tial separation p8[ 
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where the positive corrections apply to collective modes 
and the negative to relative. The collective diffusion coef- 
ficients and are enhanced by hydrodynamic cou- 
pling because fluid displaced by one sphere entrains the 
other. Relative diffusion coefficients and Djj* are sup- 
pressed, on the other hand, by the need to transport fluid 
into and out of the space between the spheres. 

Introducing a planar boundary into this system adds 
considerable complexity. The flow field around a small 
sphere located a height h above a horizontal wall is most 
easily calculated by the method of images plj ], in which 
the wall's no-slip boundary condition is satisfied by plac- 
ing a stokeslet (S), a source doublet (D), and a stokeslet 
doublet (SD) a distance h below the plane of the wall 
pj| , pO| . The flow due to this image system is described 
by the Green's function 



G^p{x ~ Rj) — -C?^g(a?— Rj) 



+ 2h*GZf> (x -Rj)- 2hG% (x - Rj ) (7) 

where Rj — ?j — 2hz is the position of sphere j's image, 
and 



G^(y) = (l-2S az )-^-( V - 
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(9) 



are Green's functions for a source dipole and a stokeslet 
doublet, respectively. The flow field set up by the image 
system (and thus by the wall's no-slip boundary condi- 
tion) entrains the sphere through b\ ai p — G^(fj — Ri) 
and decreases its mobility. Two independent modes 
emerge from this analysis, one (z) normal to the wall 
and the other (xy) parallel, with diffusivities [Oil 
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Eqs. (H) and (||) should suffice for two spheres far 
from bounding surfaces. Similarly, the spheres' mo- 
tions should decouple when the influence of a nearby 
wall dominates; Eqs. (|l0| ) and ( |ll| ) then should apply. 
At intermediate separations, however, neither set of for- 
mulas is accurate. Naively adding the drag coefficients 
due to sphere-sphere and sphere-wall interactions yields 
D' 1 (r, h) = D^ir) + D^j (h) - D^ 1 . Results of this lin- 
ear superposition approximation appear as dashed curves 
in Fig. |2[ While adequate for spheres more than 50 radii 
from the wall [(c) and (d)], linear superposition under- 
estimates the wall's influence for smaller separations [(a) 
and (b)]. 

A more complete treatment not only resolves these 
quantitative discrepancies but also reveals an additional 
surprising influence of the bounding surface on the 
spheres' dynamics: the highly symmetric and experimen- 
tally accessible modes parallel to the wall are no longer 
independent. 
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FIG. 3. Cross-sectional view of the diffusive modes for two 
spheres near a wall. Collective motion normal to the wall 
becomes increasingly coupled with relative motion parallel to 
the wall as h approaches r. Collective normal modes at large 
r cross over continuously to relative parallel modes as r de- 
creases. The dashed line at x = indicates the symmetry 
plane. 
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The combination of a neighboring sphere and two im- 
age systems contribute b\ a }(i = Gf /3 (r l - r}) + - 

Ri) + G^g(fi — Rj) to the mobility of sphere i in the a 
direction. Eigenvectors of the corresponding diffusivity 
tensor appear in Fig. ^|. The independent modes of mo- 
tion are rotated with respect to the bounding wall by an 
amount which depends strongly on both r and h. The 
experimentally measured in-plane motions clearly are 
not independent yet will satisfy Stokes-Einstein relations, 
nonetheless, with pair-diffusion coefficients D^' R (r,h) = 
Dia,ia{i",h) ± Di a ^a{r,h), where the positive sign cor- 
responds to collective motion, the negative to relative 
motion, and a indicates directions either perpendicular 
or parallel to the line connecting the spheres' centers. 
Explicitly, we obtain 

(r,h) 



Do 



1 



9 a 3 a 
16 h 4r 



1 



Do 



^ 9 a 3 a 
~ 16 h 2r 



1 - 



1 

i + 



h 4)3/2 



(i 



_ 5/2 

Ah 2 /r 2 



(12) 



(13) 



These 



up to 0(a 3 /r 3 ) and 0(a 3 /h 3 ), where £ 
results appear as solid curves in Fig. ^ 

To gauge the success of this procedure and to quan- 
tify the range over which the presence of a wall mea- 
surably influences colloidal dynamics, we computed 
the error-weighted mean-squared deviation of the pre- 
dicted diffusivities from the measured values, X^C 1 ) — 

/ [(D^\r, h) - Di,(r, hj) /6D; xpt (r, h)] * dr. Typical 
results appear in Fig. |J. The lowest-order stokeslet pre- 
sented here analysis agrees well with measurement over 
the entire experimentally accessible range. Deviations 
from the linear superposition approximation's predic- 
tions, on the other hand, are evident out to h = 15 /xm 
or 30 radii. 

The present study demonstrates that a confining sur- 
face can influence colloidal dynamics over a large range 
of separations. This influence is inherently a many-body 
effect, as demonstrated by the linear superposition ap- 
proximation's failure. Quantitative agreement between 
our measurements and a leading-order stokeslet analysis 
offers hope for future progress in understanding confine- 
ment's effects on colloidal dynamics. 
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FIG. 4. Mean-squared deviations between measured and 
predicted diffusion coefficients for relative perpendicular mo- 
tion, averaged over initial separations r. Leading-order 
stokeslet analysis yields uniformly good predictions for all ac- 
cessible values of h. Linear superposition, on the other hand, 
underestimates the wall's influence. Dashed lines are guides 
to the eye and emphasize the a/h leading-order correction to 
the linear superposition model's prediction. 
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